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Abstract 

The ferromagnetic properties of the spin-1 BEG model on finite-size Cayley tree are investigated 
using the exact recursion method. The spontaneous magnetization of the system is studied in detail 
for different values of the reduced crystal- field interaction D/J, and it is found that there is an 
unusual behavior (anti-Curie temperature) when D/J> 2.0. We also obtain the Curie temperature 
of this finite-size system. When the system size is large enough, our results will fit well with that 
in the thermodynamic limit. 
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I. INTRODUCTION 



The Cayley tree [l| which is also called Bethe lattice [2] was firstly investigated by Ku- 
rata et al. jjj. Then Domb studied the Ising model on such lattice and demonstrated that 
a Bethe- Peierls approximation is exact for the Bethe lattice {4, 5]. Over the years, the 
thermodynamic properties of the Ising system on this lattice have been extensively inves- 

n □ n ri 

igated |6|, |7|, |8|, |9j . As an expanded Ising model, the Blume-Emery-Grifliths (BEG) model 
10 . which is characterized by bilinear and biquadratic exchange interactions and crystal- 



field interaction 



12 



has played an important role in the development of the theory of 



tricritical phenomena. This model has been studied by a variety of techniques, e.g., the 

nnn 

generalized Bethe- Peierls approximation [14J, [15|, |16j, the effective-field theory 
the generalized constant-coupling approx imation theory with two parameters 



exact recursion relations method 



3, 



is 



21 



the 



22. 21 



\ Q and so on. 

All of these systems mentioned above are studied in the case of the thermodynamic 
limit. As we all know that the systems studied by the methods of experiment and numerical 
simulation are all finite. So the research on the finite-size system is much more meaningful. 
It was not until recently that a new exact expression of finite-size system for the zero-field 



magnetization was established [25]. Then the corresponding exact expression in closed form 



for the zero-field susceptibility was given by T. Stosic et al. |26|, [27|, [28(. To our knowledge 
no exact calculation has been made in the field regarding the properties of the spin-1 BEG 
model on the finite-size Cayley tree yet. 

In this paper we investigate the effect of the finite-size system on the ferromagnetic 
properties in detail. The expression of the magnetization for this system with different values 
of the reduced crystal-field interaction is derived and the Curie temperature is obtained. The 
results are compared with that of the case in the thermodynamic limit. 



II. MODEL AND FORMULATION 



At the beginning, we give a brief description of the construction of the Cayley tree. 
Starting from a single point 0, the central one of the graph [29], we add q different points 
connected to the central point which may be called "the first shell". Then each point of 
the first shell is joined to q — 1 new points. So the points of the first shell have q(q — 1) 
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nearest neighbors in total which form the second shell. The number of shells is also defined 
as generation number n while q as coordination number. If we continue in this way, the 
entire structure of the Cayley tree is formed. 

The Hamiltonian of the BEG model on the Cayley tree is given by 



(y) (ij) i * 



(1) 



where £>$(= ±1,0) is the spin at site i, the summation £\ runs over all the sites and Yluj) 
denotes summation over all the nearest-neighbor pairs. J, K, D and h describe the bilinear 
exchange, biquadratic interaction, crystal-field (or single-ion anisotropy) interaction and 
external magnetic field, respectively. This Hamiltonian was originally proposed to explain 
the phase separation and superfluidity in 3 He- 4 He mixtures 101 ] . 
The partition function of the above system can be written as 



Z = ^exp (-(3H) 



J^exp 



P [ JY,S 1 S 3 + K^2S*S* - Dj^S* + hJ2 Si 

(ij) (ij) i i 



(2) 



where (3 = 1/ksT, ks is Boltzmann constant and T is the absolute temperature. The 
summation ^ 5 in Eq. §2§ goes over all spin configurations of the system. 

Without loss of generality, we consider a Cayley tree of n generations with branch number 
B = 5 (coordination number minus one). Then, the n-generation branch consists of N n = 
5 " + 4 1 ~ 1 spins, while the 0-generation being a single spin. Let zi^ and Zn^ be the 

partial partition functions of the system, with the central spin takes values +1, and -1 
respectively. Based on Eq. (j2j), we can obtain the recursion relations |6|] as 



7 (±) 

J n+1 



and 



7 (o) 

J n+1 



( Z (+) + Z (o) + z (-)) 5 . 

\ n 1 n 1 n / 



(3) 



(4) 



The partition function Eqs. (j3J) and (jlj) can be differentiated with respect to field, thus 
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the recursion relations for the field derivatives of the partition function are easily written as 

F>7 {±) 

^±1 = ±e -/^ ( e ±PJ + PK z t+) + Z (0) + e TW Z H) 5 
CJ L) f v 

/ A7(+) Wo) W-)' 

I d/3h d(3h d(3h 

Wo) . /W+) Wo) W-)' 

S/9/i " 5( n + " + " j ^/3/i + d(3h + d(3h 

Starting from a single spin (0-th generation branch) we have 



Zj — e , Zjq 



and 



-/3D±0h <9Zq^ _ 

Then the magnetization of a site in the Cayley tree can be written as 

<m> - = «^ W (6) 

The exact recursion relations and the magnetization expression allow us to study the 
thermodynamic behavior of the system in detail. In next section, we will give the numerical 
results of the magnetization. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

In the following, we investigate numerically the magnetization of the BEG model on this 
finite-size Cayley tree. Based on the above Eqs. ([IH6]), taking the limit h — > 0, we make a 
detailed calculation to the magnetization for various strengths of the interaction J, K and 
D and for several system sizes n. 

The magnetization as a function of temperature is shown in Fig. 1 for various reduced 
values of crystal-field interaction D/J while K/J = 1, n = 12. The curve corresponding to 
D/J = 2.0 separates the curves into two different trend. In the case of D/J > 2.0, there 
is an anti-Curie temperature T aC [yj], i.e., a temperature below which the magnetization 
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vanishes, in addition to the Curie temperature Tc- As D/J decreases, the anti-Curie tem- 
perature drops to lower values and the maximum magnetization becomes larger. We can see 
that the anti-Curie temperature vanishes in the case of D/J = 2.0 and the magnetization 
decreases steadily from its saturation value 0.58 to zero with growing temperature. In the 
case of D/J < 2.0, the maximum magnetization of curves increases and it approaches to 
unity at zero temperature for D/J = 1.5. Moreover the magnetization curve would keep 
the conventional shape for even smaller D/J values that is not shown in this figure. An 
approximative value of the deduced Curie temperature T c = 1.4 can also be obtained for 
this two-order transition. 

In Fig. 2 we present the magnetization as a function of temperature for several system 
sizes n = 3,6,9, 12 and different values of D/J. It is seen that the curves exhibit a slow 
decay of magnetic ordering with the increase of the system sites. The larger the system size, 
the sharper the magnetization curve is. If the system sizes could get large enough value, 
the magnetization curve of finite-size system would become close to that in thermodynamic 
limit. The unusual shape of the figure is due to the complicated interactions including the 
bilinear exchange and biquadratic exchange. This result is in agreement with the system in 



thermodynamic limit which is studied by K. G. Chakraborty et al. 
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3l|. 



IV. CONCLUSIONS 



In this paper, using the exact recursion relations, we present an exact calculation for 
the spontaneous magnetization of the spin-1 BEG system on the finite-size Cayley tree. 
The magnetization properties are studied in detail for different values of the crystal-field 
interaction and system size of the Cayley tree. It is shown that the magnetization curves 
exhibit some unusual features including an anti-Curie temperature as the variance of the 
strength of the crystal-field interaction. The Curie temperature is also obtained. The curves 
would become close to that in the thermodynamic limit with the increase of the system size. 
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FIG. 1: The magnetization M as a function of T for K/J = 1, n = 12 and 5 = 5. The 
successive curves from (a) to (b) are for D/ J = 2.3, 2.1, 2.0, 1.9, 1.5 and 1 respectively. There 
exhibits a two-order phase transition. 

FIG. 2: The magnetization versus temperature T for different D/J and n = 3,6.9,12. 
The successive curves from (a) to (f) are for D/J = 2.3,2.1,2.0, 1.9, 1.5 and 1 respectively. 
The slope of each subgraph curve steepens with the increase of the system sizes n. 



8 




FIG. 1: The magnetization M as a function of T for K/J = 1, n = 12 and B = 5. The successive 
curves from (a) to (b) are for D/J = 2.3,2.1,2.0,1.9,1.5 and 1 respectively. There exhibits a 
two-order phase transition. 
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FIG. 2: The magnetization versus temperature T for different D/J and n = 3,6.9,12. The 
successive curves from (a) to (f) are for D/J = 2.3,2.1,2.0,1.9,1.5 and 1 respectively. The slope 
of each subgraph curve steepens with the increase of the system sizes n. 
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